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The  finite  basic  solution  method  is 
employed  in  this  paper  for  calculating  the 
subsonic  lift  and  drag  characteristics  of  a 
back-swept  wing  with  wing-tip  mounted 
winglets.  Calculations  have  also  been 
performed  on  its  pitching  moment  and  wing 
root  bending  moment.  Through  calculations  of 
various  wing-tip  mounted  winglet  configura¬ 
tions,  the  general  rules  of  wing-tip  mounted 
winglet  aerodynamic  characteristics  have  been 
analyzed.  A  comparison  between  calculated 
results  and  experimental  data  reveals,  that 
this  method  can  satisfy  the  needs  for  the 
preliminary  design  and  configuration 
selection  of  a  wing-tip  mounted  winglet. 
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I .  Preface 

A  few  papers  concerning  theoretical  calculations  of  wing-tip 
mounted  winglets  have  been  seen  in  foreign  journals  in  recent 
years.  Yet,  basically  no  experimental  verifications  have  been 
conducted.  Therefore,  the  extent  of  their  practical  applications 
cannot  be  determined.  At  present,  in  China  and  abroad,  research  of 
configuration  selections  for  wing-tip  mounted  winglets  generally 
employs  experimental  methods.  This  is  both  expensive  and  time 
consuming.  Moreover,  only  a  few  types  of  winglet  are  used  in 
experiments  and  their  configurations  are  relatively  simple.  Based 
on  the  existing  data,  the  general  rules  of  wing-tip  mounted 
winglet  aerodynamic  characteristics  are  still  inconclusive.  This 
paper  presents  a  theoretical  calculation  method  utilizing  the 
finite  basic  solution  method  to  calculate  the  subsonic  lift  and 
drag  characteristics  of  a  back-swept  wing  with  wing-tip  mounted 
winglets.  Calculations  have  also  been  performed  for  the  pitching 
moment  and  the  wing  root  bending  moment.  Thereby  the  general  rules 


of  the  wing-tip  mounted  winglet,  geometrical  dimensions,  and  the 
winglet  rake  angle,  etc.  can  be  studied.  In  order  to  verify  the 
applicabilities  of  this  calculation  method,  force  and  pressure 
experiments  were  conducted  on  some  of  the  calculated  examples.  The 
calculated  results  basically  conform  well  with  the  experimental 
data;  therefore,  this  calculating  method  has  practical  signifi¬ 
cance  in  the  preliminary  design  and  configuration  selection  of  the 
wing-tip  mounted  winglets. 

II.  Calculation  Method  and  Results 

In  the  finite  basic  solution  method  employed  in  the  paper, 
the  equal  deflection  angle  method  is  used  to  plot  the  spanwise 
grid  and  to  determine  the  locations  of  the  spanwise  control 
points;  for  chordwise,  the  equal  partial  method  determines  the 
grid.  Grids  are  laid  out  on  the  wing  and  wing-tip  mounted 
winglets.  A  compressible  horseshoe  vortex  is  laid  out  on  each  grid 
and  the  attached  vortex  is  placed  at  the  1/4  chord  line  in  the 
grid.  Two  free  vortices  are  extended  along  the  streamline  to 
infinity.  Chordwise  locations  of  the  control  points  are  set  at  the 
3/4  chord  line  of  the  grid.  The  flow  boundary  conditions  at  each 
control  point  of  the  main  wing  and  wing-tip  mounted  winglets  are 
satisfied.  A  combined  flow  field  method  is  employed  for 
calculating  induced  drag  Cl,  21. 

The  two-dimensional  sizes  of  back-swept  main  wing  and  winglet 
calculated  in  this  paper  are  shown  in  Fig.  1. 


Fig.  1.  The  two-dimensional  sizes  of  swept 
main  wing  and  winglet. 

Key:  (l)  The  main  wing;  (2)  The  winglet. 

During  the  calculation,  wing  and  winglet  are  divided  into  5 
columns  in  the  chordwise  direction;  and  into  a  total  of  10  rows  in 
the  spanwise  direction  with  7  rows  occupied  by  the  main  wing  and  5 
rows  by  the  winglet.  If  the  winglet  is  cut  in  the  chordwise  direc¬ 
tion.  it  is  divided  into  five  strips.  Each  strip  may  possess  a 
difrerent  angle  of  rake  so  that  various  different  winglet  arrange¬ 
ments  can  be  formed.  A  total  of  14  winglet  arrangements  have  been 
calculated  in  the  paper.  The  14  winglet  arrangements  and  their 
three-dimensional  schematic  diagrams  are  shown  in  Fig.  2. 
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Fig.  2.  The  three-dimensional  figures 
or  winglet  arrangements. 

The  degrees  of  rake  angle  of  the  5-strip  winglets  are  given 
in  the  parentheses  in  Fig.  2.  The  order  of  the  rake  angles  is  from 
the  leading  to  the  trailing  edges.  The  wing-chord  plane  is  the 
datum  plane  for  angle  measuring.  Upward  rake  is  positive  and  down¬ 
ward  rake  is  negative.  Arrangement  1  (X,  X.  X.  X.  X)  shown  in 
Fig.  2  is  a  main  wing  without  mounted  winglet.  In  order  to 
facilitate  the  comparison  of  aerodynamic  characteristics  for  each 
winglet  arrangement,  the  area,  the  span  to  chord  ratio,  the  mean 
aerodynamic  chord  and  the  span  length  of  the  main  wing  are  used  as 
the  reference  area  and  reference  lengths  during  calculations  of 
aerodynamic  coefficients  for  various  arrangements.  The  aerodynamic 
characteristics  for  14  winglet  arrangements  obtained  from  calcula¬ 


tions  are  listed  in  Table  1. 


Table  1 . 


The  result  ot  calculation 
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Key:  (1)  Arrangement. 


Comparing  the  data  of  various  winglet  arrangements,  the 
following  observations  can  be  made! 


1)  No  matter  what  kind  of  winglet  is  added  to  the  main  wing. 


its  C*  value 
a  r  rangement 
i mproved  as 
and  are 


is  always  somewhat  higher  than  that  of  main  wing 
1.  The  induced  drag  characteristics  are  somewhat 
the  5  value  is  decreased.  However,  values  of  its 
also  somewhat  increased  as  value  of  C"  is  increased. 


2)  It  can  be  seen  from  the  trends  of  arrangements  2^5  that, 
the  smaller  the  rake  angle  of  the  one-piece  winglet  is,  the  larger 
C?  is  and  the  smaller  S  is.  The  aerodynamics  coefficients  for 
arrangements  3  and  4  are  identical,  which  shows  that,  regardless 
of  whether  the  winglet  rake  is  positive  or  negative,  the  aero¬ 
dynamic  characteristics  of  a  one-piece  winglet  are  basically 
identical  so  long  as  the  degrees  of  rake  angle  remain  the  same. 
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3)  For  the  two-piece  winglet  arrangement  6,  the  aerodynamic 
performances  of  C“  and  <S  are  not  as  good  as  that  of  the  wing 
elongating  winglet  arrangement  5.  However,  it  is  better  than  the 
one-piece  winglet  arrangements  2~4. 

4)  For  a  mu  1 t i p 1 e -p i ece  winglet  arrangement  (rake  angles  of 
the  5  winglet  strips  are  all  different),  its  C“  and  ds  are  better 
than  other  winglet  arrangements.  The  values  of  C“  and  8  are 
similar  for  the  various  mu  1 t i p 1 e - w i ng 1 e t  arrangements.  Therefore, 
from  the  standpoint  of  increasing  the  lift  and  decreasing  the 
induced  drag,  a  mu  1 t i p 1 e- w i n g 1 e t  arrangement  is  a  better  selected 
ar rangement . 

The  above  four  concluding  remarks  are  consistent  with 
conclusions  for  rectangular  wing  in  reference  [1],  which 
illustrates  that,  regardless  of  a  rectangular  or  a  back-swept 
wing,  the  trends  of  aerodynamic  characteristics  changes  are 
identical  after  a  wing  tip  winglet  is  mounted  on  it. 

111.  Experimental  Method  and  Results 

The  experiments  described  in  this  paper  were  performed  in  a 
NH-2  low-speed  wind  tunnel  at  Nanjing  Institute  of  Aeronautics  and 
Astronautics.  The  wind  speed  of  experiments  was  40  m/s.  The  main 
wing  model  used  was  a  wooden  semi-mould  model.  Holes  for  pressure 
tappings  were  drilled  on  the  wooden  main  wing.  The  main  wing  wing- 
tip  was  mounted  with  winglets.  Winglets  were  made  of  aluminum 
plate  of  1.5  mm  in  thickness.  No  pressure  tapping  holes  were 
drilled  on  the  winglets.  Thus,  for  that  model,  measurements  of 
pressure  distributions  were  made  only  on  the  main  wing.  In  order 
to  save  both  experimental  time  and  expenditure,  only  winglet 
arrangements  1~9  were  investigated.  Experimental  work  included 
force  tests  and  some  pressure  tests.  The  angle  of  attack  for  aii 
experiments  ranged  from  4° ~  24° .  The  main  characteristics  data 
obtained  from  force  test  experiments  are  listed  in  Table  2. 
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Tab  i  e  2 . 


The  result  of  experiment 
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Key:  i  1  )  Arrangement. 


It  can  be  seen  from  the  comparison  between  data  listed  in 
Tables  1  and  2  that  the  magnitudes  of  the  calculated  values  for 
C? .  mf*  ,  and  m$»  are  basically  similar  to  those  of  the  experi¬ 
mental  data;  except  that  the  calculated  values  of  <£  are  markedly 
smaller  than  those  of  the  experimental  data.  Experimental  findings 
and  calculated  results  of  the  aerodynamic  characteristics  for 
various  winglet  arrangements  are  found  to  be  consistent.  In  other 
words,  the  aforementioned  four  concluding  remarks  which  were 
derived  from  calculations  have  all  been  verified  in  the  experi¬ 
ments,  except  for  the  following  two  minor  discrepancies  between 
the  experimental  values  and  the  calculated  results. 

1)  The  experimental  value  of  £  in  winglet  arrangement  4  is 
less  than  that  of  winglet  arrangement  3. 

2)  The  experimental  value  of  S  in  winglet  arrangement  6  is 
less  than  that  of  winglet  arrangement  5. 

Calculations  and  experiments  have  all  proved  that  the  values 
of  £  for  mu  1 t i p 1 e - p i ece  winglet  arrangements  7~9  are  less  than 
those  for  2  ~ 6  winglet  arrangements.  This  indicates  that  the 
induced  drag  characteristics  of  mu  1  t i p 1 e - p i ece  winglets  are  better 


Equations  (.19)  — '  (21)  are  further  expressed  as  calculation 
r  orms .  We  have 
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At  any  given  time  instant,  if  values  of  dP/dt  ,  dv/dt  •  and  dp/di 
in  equations  (22)/^/(24)  are  known,  the  above  set  of  equations  can 
be  regarded  as  an  ordinary  differential  equation  set  with  respect 
to  x.  The  solution  can  now  be  obtained  by  using  the  Triple-Loop- 
Iteration  Algorithm  described  above.  After  the  values  of  p.  v,  and 
p  at  all  node  points  in  the  axial  direction  are  obtained  for  a 
certain  time,  instant,  they  are  substracted  by  their  corresponding 
values  at  the  previous  time  instant  and  then  divided  by' the  time 
interval  between  the  two  time  instants,  thereby  the  mean  change 
rates  for  the  P,  v,  and  p  values  with  respect  to  time  at  any  node 
point  between  those  two  time  instants  are  calculated.  With  the 
mean  change  rates  as  the  approximation  for  dp/dt  .  dv/dt .  and  dp/dt 
and  substituting  them  into  equations  !22)~(24),  these  equations 
can  then  be  solved  by  repeating  the  aforementioned  method.  Thus, 
by  repeatedly  alternating  the  calculations  with  respect  to  time 
and  space  independent  variables,  the  dp/dt,  dv/dt,  and  dp/dt  values 
are  corrected  through  iteration.  After  a  certain  accuracy  order  is 
achieved  the  calculations  for  the  next  time  instant  are  initiated 
and  the  whole  process  is  continued  until  the  rocket  motor 
operation  is  completed. 


By  rewriting  equations  <.  2  2  )  ~  t  2  4  >  in  a  vector  form,  we  have 
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The  essence  of  Time-Space  Alternate  Algorithm  is 
a  set  of  partial  differential  equations  into  a  set  or  ordinary 
differential  equations.  With  respect  to  time  or  space  independent 
variables,  a  higher  order  of  accuracy  for  the  solution  is  obtained 
by  using  this  algorithm.  In  reference  [6],  Burstein  employed  a 
similar  algorithm  for  solving  hyperbolic-type  equations  and 
obtained  a  third  order  accuracy  for  both  time  and  space  indepen¬ 
dent  variables.  The  Time-Space  Alternate  Algorithm  employed  in 
this  paper  not  only  achieves  a  fourth  order  accuracy  with  respect 
to  space  independent  variable  but  also  provides  a  good  link 
between  solution  methods  for  the  unsteady  equations  and  the 
quasi -steady  equations. 

In  the  Time-Space  Alternate  Algorithm.  P,  v.  and  p  are 
selected  as  dependent  variables.  First,  equations  (  1  ) '■'•'(  3  )  are 
rearranged  as 
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parameters  at  time  t„  should  first  be  obtained  from  grain 
calculation.  Within  time  interval  A  t .  the  thickness  Ay,  of  the 
burned-out  grains  is  calculated  by  the  following  formula: 

A  y/-»r,A* 

where  r,  is  the  average  combustion  rate  within  the  time  interval 
At.  This  value  is  dependent  upon  the  combustion  rate  at  time  t« . 
However,  the  combustion  rate  is  also  affected  by  all  combustion 
gas  parameters  and  consequently,  it  is  also  affected  by  Ay, . 
Therefore,  it  must  be  solved  for  through  the  successive  iteration 
method.  We  have  employed  a  single-prediction  and  multiple- 
correction  approach.  Details  of  the  iteration  formulae  are: 

(w-i,  2,  •••) 

rr  a.)  -  0,0-.,) + rrv  cox  o -i,  2,  •••) 

When  employing  the  Merson  Method  to  conduct  numerical 
integration,  it  requires  grain  geometrical  parameters  at  certain 
points  between  two  node  points  and  also  their  derivative  values 
with  respect  to  X.  These  data  are  obtained  by  the  interpolation  of 
LaGrange  quadratic  polynomial  equation  and  by  the  corresponding 
numerical  differential  formulae. 

3.  Time-Space  Alternate  Algorithm 

The  governing  equations  ( 1  )  <~'  (  4  )  of  a  one-dimensional 
unsteady  internal  ballistic  performance  analysis  are  the  first 
order,  linear,  hyperbolic  partial  differential  equations.  Due  to 
the  limitations  in  the  difference  method's  stability  condition, 
the  time  increment  must  be  set  very  small  if  the  characteristic 
line  method  is  used  to  solve  these  equations.  For  a  typical  rocket 
motor,  ten  thousand  or  even  more  iterations  are  required.  This  is 
verv  inconvenient  and  sometimes  even  unacceptable. 
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After  determining  boundary  conditions  of  pa  and  mG  at  the 
motor  head,  the  Merson  Method  is  employed  to  calculate  the  initial 
condition  problem  of  equations  (10)  and  (11).  Then,  a  correction 
equation  is  applied  to  correct  the  pc  value  in  accordance  with  the 
exit  boundary  conditions.  For  this  purpose,  the  generation  rate  of 
the  combustion  gas  m„  and  the  combustion  gas  exhaust  rate  m,  are 
given  by 

Bti*  mi  +  itij 

frit  j^i 

where  ma  -  the  generation  rate  of  the  combustion  gas  on  the 
tail-end  surface  of  the  grain.  The  flow  from  the  cross-section  at 
grain  exit  to  the  nozzle  throat  is  also  considered  quas i - s teady . 
Consequent  1 y ,  m»  should  be  equal  to  m„  .  After  one  iteration,  if 
the  calculated  values  of  mt  and  m„  are  not  equal,  it  indicates 
improper  p„  guess  value  and  the  value  has  to  be  corrected.  The 
first  correction  for  p0  is  performed  by  the  following  'proportions  1 
formula: 


After  the  number  of  iteration  exceeds  two,  an  interpolation 
formula  can  be  used  for  p„  correction. 

Let  — m*.  thus. 


„(»•!)  _ 

‘pru+ Am1-1’  Am0.-t)_AV-<r 

(  n  -3,  4,  — ) 


teration  is  teminated  when  Am  is  less  than  a  specified  value. 


When  proceeding  with  the  solution  at  time  t.  after  the 


solution  at  time  t 


has  been  obtained,  the  geometrical 


The  solution  method  for  the  set  of  internal  ballistic 
equations  is  a  kind  of  boundary  condition  problem  and  is  not  an 
initial  condition  problem.  Thus,  an  initialization  guess  method  i 
employed.  First,  assuming  the  guess  value  of  pressure  pD  at  the 
motor  head  and  at  time  ta  can  be  estimated  in  accordance  with  the 
following  formulae: 


\  Pronto,  2  (  Y  +  1  )  1  }  -frr 

l  A,  y  I  CXi)  ) 


where 


4(l  -X») 


(i+Xl)*"(l -y+-fxiJ£(X) 


The  above  \  ,  is  determined  by 


<  12  > 


<  1 


<  1  a  ) 


At  time  t, ,  the  initial  guess  value  of  p0  is  taken  from  the  value 
of  p0  which  satisfies  continuity  at  time  t„ .  At  time  t2  and  other 
later  time  instants,  initial  guess  values  of  pa  are  obtained  bv 
the  following  proportion  using  the  p„  values  at  two  previous  time 
i nstants ;  i . e . 


/%(*-)-/>. (™-2.  3’  -) 
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At  any  time  instant,  the  mass  flow  rate  across 
section  0  should  be  equal  to  the  generation  rate  of 
on  the  combustion  surface  at  the  motor  head.  Hence, 


the  cross- 
combus  t i on 


m^AttPraxpl 


>.  16  ) 


If  a  ^  a ,  P  <£p,,  A  Ap<£Ar  (within  the  residence  time  of  the 
combustion  gas) ,  the  operation  process  of  the  rocket  motor  can  be 
considered  quas i -steady.  The  above  equations  can  be  simplified  as 


_«L(p  ^)  =  P,r-^  (5) 

-j-(Pv>Af+pA')  (6) 

^[pvAf(cr+H))=Prrd^Hr  (7) 

p =PRT  (8) 

2.  Tr i p 1 e-Loop- I terat i on  Algorithm 

With  respect  to  a  one -d i mens i ona  1  quas i - s teady .  adiabatic 
flow,  energy  equation  (7)  in  the  above  governing  equations 
(5)~(8)  can  be  simplified  as 

crT+x>'ll-crT.  (9) 

Consequently,  only  the  continuity  and  momentum  equations 
among  the  governing  equations  are  expressed  in  differentiated 
form.  If  the  mass  flow  rate  m  and  the  pressure  p  are  selected  as 
the  calculating  variables,  the  continuity  and  momentum  equations 
can  be  written  as  follows: 


dm 

dx 


=  PfT 


dA* 

dx 

1 


dp  _ 

dx  Ari^-v1)  (h 


---  v  C2o*+(  Y  -  1  | 


(  10 ) 

(11) 


By  using  the  energy  equation  and  equation  of  state  and  other 
well-known  equations,  all  the  unknown  parameters  can  be  expressed 
as  functions  of  m  and  p. 


internal  ballistic  calculations  are  closely  related  and  must  not 
be  separated,  a  general  two-dimensional  grain  calculation  method 
introduced  by  the  author  is  employed. 


II.  Tr i p 1 e-Loop- I t e r a t i on  and  Time-Space  Alternate  Algorithm 


1.  Rocket  motor  model  and  its  one - d i mens i ona 1  internal  ballistic 


governing  equations 


Figure  1  illustrates  a  selected  solid  propellant  rocket 


motor . 


i 


Fig.  1.  Schematic  of  solid  propellant  rocket  motor 


For  the  simulation  of  a  combustion  chamber  flow  field,  the 
following  assumptions  have  been  made:  The  combustion  gas  is  an 
ideal  gas;  the  fluid  flow  through  the  hollow  section  of  a  rocket 
motor  cylinder  is  a  one -d i mens i ona 1  flow;  the  fluid  flow  process 
is  a  frictionless,  adiabatic,  non- react i ve,  homogeneous  and  gas- 
phase  hyd rodynam ic  process.  The  projection  of  the  momentum  of 
combustion  gas,  which  is  being  generated  on  the  combustion 
surface,  on  the  x-axis  is  also  assumed  to  be  zero.  The  governing 
equations  of  the  fluid  flow  in  the  rocket  motor  can  be  derived  as 

~jr(PAr) +~h~(  'M*)“p,r^sr  a) 

-fr (pvAf)  +  ~hr  -  p  ~  ( 2 ) 

- Hpa '  (rr + £)] + [pvA>  Or- + ")] = p'r  H>  ( 3  > 

P-PRT 


The  main  purpose  for  predicting  the  solid  propellant  rocket 
motor's  internal  ballistic  performance  is:  based  on  the  original 
data  of  rocket  motor  structures,  grain  geometrical  shapes  and 
dimensions,  propellant  characteristics,  etc.  to  calculate  the 
distributions,  along  the  length  of  the  combustion  chamber,  of 
parameters  such  as  combustion  gas  pressure  p,  flow  velocity  v. 
density  P,  temperature  T,  etc.  during  the  motor  operation  period 
and  their  changes  with  respect  to  time  t. 

Throughout  the  years,  there  has  been  great  progress  in  the 
development  of  internal  ballistic  calculation  algorithms.  The 
earlier  internal  ballistic  calculation  algorithms  were  established 
on  the  simplified  bases  of  zero-dimensional  and  quasi-steady  flows 
C13.  Later,  one-dimensional  flow  model  was  employed  as  the 
volumetric  filling  coefficient  became  larger  and  the  erosion 
combustion  phenomenon  became  very  distinctive.  Based  on  the  one¬ 
dimensional  fluid  flaw  model  analysis  and  further  simplifications 
to  the  quasi-steady  flow,  a  set  of  ordinary  differential  equations 
is  obtained  C2~5].  Thus,  the  Runge-Kutta  Method  or  Merson  Method 
can  be  applied  for  solutions  C7],  If  the  flow  is  one -d i mens i ona 1 
and  unsteady,  the  governing  equations  for  the  rocket  motor  flow 
calculations  become  a  set  of  first  order,  linear,  hyperbolic 
partial  differential  equations.  Under  these  circumstances,  a 
solution  can  be  obtained  by  the  finite-difference  and  charac¬ 
teristic  line  methods  C8]. 

Emphasis  of  this  paper  is  placed  on  the  introduction  of  a 
calculation  algorithm  generally  applicable  to  both  quasi-steady 
equations  and  unsteady  equations;  i.e.  the  Time-Space  Alternate 
Algorithm.  The  algorithm  was  developed  on  the  basis  of  a  Triple- 
Loop-  1 teration  Algorithm.  Since  the  grain  calculations  and  the 
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a,  Burning  speed  coefficient 

o,  Critical  speed  of  sound 
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Two  algorithms  for  prediction  of  solid 
propellant  rocket  motor  performance  are  intro¬ 
duced,  and  the  emphasis  is  placed  on  the  Time- 
Space  Alternate  Algorithm.  A  two-dimensional 
grain  general  calculation  method  introduced 
by  the  authors  for  the  calculation  of  grain 
is  employed.  Thus,  the  method  presented  in 
this  paper  for  predicting  the  internal 
ballistic  performance  has  the  generality  of 
applicabilities.  A  comprehensive  computer 
program  for  the  above  method  has  been  written 
and  applied  to  calculating  the  performance  of 
three  different  solid  propellant  rocket 
motors.  The  calculated  results  are  consistent 
with  the  experimental  data. 


It  can  be  observed  from  the  comparison  between  the  theoreti¬ 
cal  calculations  and  experimental  data  that  said  theoretical  cal¬ 
culating  method  can  be  employed  to  conduct  preliminary  evaluation 
of  various  winglet  arrangements,  and  it  has  practical  significance 
in  the  preliminary  design  and  configuration  selection  of  the 
wing-tip  mounted  winglets. 
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Fig.  4.  The  experimental  curve  of  pressure  distribution. 

IV.  Concluding  Remarks 

It  is  concluded  from  the  above  analyses  that,  the  smaller  the 
rake  angle  of  the  one-piece  winglet  is.  the  larger  the  Cy  Df  ^e 
main  wing  is  and  the  smaller  &  is.  However,  the  length  of  wing 
span  and  the  bending  moment  at  the  wing  root  are  also  increased. 
Therefore,  in  foreign  countries,  winglets  with  rake  angle  of  about 
75°  are  often  selected  when  an  existing  transport  requires  winglet 
mounting. 

The  shape  of  a  mu  1 t i p 1 e - p i ece  winglet  resembles  that  of  an 
eagle's  wing-tip.  It  has  the  best  induced  drag  characteristics. 
However,  the  utilization  of  this  type  of  winglet  arrangement  on 
the  aircraft  often  creates  certain  structural  difficulties.  There¬ 
fore,  two-piece  winglets  are  usually  mounted  on  transport  aircraft 
in  foreign  countries. 

Of  course,  the  selection  of  a  better  arrangement  for  mounting 
wing-tip  winglets  should  be  based  on  specific  circumstances  and 
various  factors;  it  can  only  be  determined  after  a  comprehensive 
analysis  of  its  advantages  and  disadvantages. 


The  advantage  of  lift  coefficient  for  several  primary 
winglets  and  their  relationships  can  be  seen  in  the  experimental 
curves  in  Fig.  3  after  winglets  are  mounted  on  a  main  wing.  It 
should  be  pointed  out  that  the  relationships  shown  in  Fig.  3  for 
lift  coefficients  of  various  winglet  arrangements  no  longer  apply 
after  the  angle  of  attack  exceeds  8°.  The  characteristics  of  lift 
to  drag  ratio  also  display  a  similar  situation. 

The  calculations  in  this  paper  and  reference  [1]  indicate 
that,  after  winglets  are  mounted,  an  increase  in  circulation  flux 
at  the  main  wing  outboard  region  is  the  main  reason  for  Cy  to 
increase.  The  experiments  of  pressure  distribution  shown  in  this 
paper  verify  that,  after  the  winglets  are  mounted,  the  pressure 
distribution  on  the  inboard  region  of  a  main  wing  is  basically  the 
same  as  that  on  one  without  winglets.  The  cross-sectional  pressure 
distributions  at  86.5%  main  wing  span  are  measured  after  several 
primary  winglets  are  mounted  and  the  experimental  curves  are  shown 
i n  Fi g .  4. 

The  above  experimental  curves  of  pressure  distributions  show 
that  the  pressure  distribution  on  the  outboard  region  of  the  main 
wing  is  markedly  increased.  This  indicates  that  the  conclusions  of 
theoretical  calculations  are  correct. 
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than  those  of  the  one-piece  and  two-piece  winglets.  However. 


values  of  K. .  .  for  winglet  arrangements  7~ 9  are  often  less  than 
those  for  winglet  arrangements  2~'6.  From  the  viewpoint  or 
obtaining  a  larger  maximum  lift  to  drag  ratio,  the  multiple-piece 
winglet  is  not  necessarily  a  better  arrangement.  The  reason  is 
that  the  C„0  values  for  mu  1 t i p 1 e- p i ece  winglets  are  larger  (see 
Table  2)  and  that,  even  though  their  induced  drag  characteristics 
are  better,  the  larger  C,o  values  still  cause  the  «...  value  to 
decrease . 

Figure  3  illustrates  several  experimental  curves  or  lift 
coe  f  f i c i en t . 


Fig.  3.  The  experimental  curve  of  lift  coefficient. 
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In  order  to  facilitate  the  discussion  of  the  stability 
problem  of  the  difference  method,  the  above  equation  is  generally 
simplified  as: 

*2  A  d 2 

*x  ~jT  (26) 

Here,  Aa  is  a  constant  matrix. 

When  the  Time-Space  Alternate  Algorithm  is  applied  for  the 
solution  of  model  equation  (26),  it  is  equivalent  to  applying  the 
difference  method. 


_  A  [K'~ 

~  6  l  T 


-+  4 


2**} 
**/♦  i/i 


l/t 


The  finite-difference  method  is.  with  respect  to  the  model 
equations,  an  unconditionally  stable  scheme.  It  has  been  verified 
through  actual  calculations  that  the  above  conclusion  is  also 
applicable  to  the  actual  equations. 


The  Tr i p 1 e-Loop- I terat i on  Algorithm  is  employed  in  con¬ 
sideration  of  the  solutions  with  respect  to  the  space  variable  x 
in  the  Time-Space  Alternate  Algorithm.  Thus,  only  the  difference 
between  the  present  method  and  the  quasi-steady  equations  solution 
method  is  introduced  here. 


The  initial  conditions  required  for  the  Time-Space  Alternate 
Algorithm  include  initial  values  of  P,  v,  and  p  at  all  node  points 
as  well  as  the  initial  values  of  iP/M,  dv/dt  ,  and  dp/M  .  When  the 
calculations  start  immediately  after  the  rocket  motor  ignition 
process,  the  above  initial  conditions  can  be  determined  by  the 
results  calculated  from  the  ignition  process.  If,  however,  the 
calculations  start  directly  from  the  steady  state  period  of  the 
rocket  motor,  we  may  assume  that  the  flow  field  inside  the  com¬ 
bustion  chamber  reaches  steady  state  instantaneously.  Thus,  the 
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values  of  aP/at  ,  av/at ,  and  ap/at  at  all  node  points  can  be  initial  1 
set  as  zero  and  the  values  of  P.  v.  and  p  at  al I  node  points 
obtained  by  the  Triple-Loop-Iteration  Algorithm  are  used  as 
initial  values  in  the  calculations. 

Among  the  boundary  conditions  at  the  rocket  motor  head,  the 
method  for  determining  the  po  value  is  the  same  as  the  Triple- 
Loop-Iteration  Algorithm.  The  values  of  and  v0  ,  however,  are 
determined  by  using  the  following  formulae: 

Pt-pt/RT> 

_  - 

*  P.^. 

At  time  t«  ,  the  values  of  AP/A t.  Aw/ A/,  and  Ap/Ai  between  t,  .  2 
and  t.  .  ,  are  used  as  the  guess  values  for  aP/ai  ,  av/at,  and  ap/at  at 
al  1  node  points.  Having  calculated  values  of  p.  v,  and  p  at  all 
node  points  for  time  t«  ,  their  corrections  are  subsequent  1 y . made . 
The  correction  formula  is 

f  aa  V1"’  _  ara.)-g,ft-,) 

\  Jl  f.-f-, 

(w-1,  2,  • »  n  —l,  2,  •••) 

Corrections  can  be  made  repeatedly  until  the  difference 
between  the  values  of  P's,  v's,  and  p's  obtained  from  two 
successive  iterations  is  less  than  a  specified  value. 


III.  Calculation  of  Grain 

The  purpose  of  the  grain  calculation  is  to  obtain  the  grain 
geometrical  parameters  required  for  the  instantaneous  internal 
ballistic  calculations.  For  different  types  of  grain  calculation, 
not  only  the  forms  of  the  required  grain  geometrical  parameters 
are  different,  the  correlation  between  grain  and  internal  ballis¬ 
tic  calculations  may  also  be  different.  For  a  zero-dimensional 
internal  ballistic  calculation,  the  required  grain  geometrical 
parameters  generally  correspond  with  those  values  of  the  grain 
combustion  area  at  a  given  burned-out  thickness.  Given  this,  the 
grain  and  the  internal  ballistic  calculations  can  be  performed 
separately.  On  the  other  hand,  for  the  internal  ballistic 
calculation  of  a  one -d i mens i ona 1 ,  variable  cross-sectional  tunnel, 
the  required  grain  geometrical  parameters  generally  correspond 
with  those  values  of  combustion  circumference  and  tunnel  area  tor 
a  partially  burned-out  thickness  at  all  node  points.  In  view  or 
the  existence  of  erosion  burning,  the  grain  and  the  internal 
ballistic  calculations  are  closely  related  and  cannot  be 
separated.  Thus,  calculation  must  be  performed  iteratively  and 
a  1 terna  te 1 y . 

The  method  introduced  in  this  paper  for  the  internal  ballis¬ 
tic  calculation  is  coordinated  with  the  calculation  of  the  grain. 

A  general  two-dimensional  grain  calculation  method  [9]  is 
presented  by  the  authors.  Based  on  this  method,  the  solution 
scheme  can  be  extended  to  solving,  for  example,  three-dimensional 
calculation  problems  of  the  double-arc,  transitionally  closed 
motor  head  and  the  unrestrictive  combustion  tail  of  the  tail-end 
surface  located  on  the  oval  portion  following  the  motor. 
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IV.  Sample  Calculations  and  Discussion 


A  comprehensive  computer  program  for  the  above  mentioned 
methods  was  written  and  applied  to  the  calculation  of  the  perform¬ 
ance  of  three  different  real  rocket  motors.  For  each  motor,  both 
the  Tr i p 1 e-Loop- 1 terat i on  Algorithm  and  the  Time-Space  Alternate 
Algorithm  were  employed.  In  examples  1  and  3,  since  the  internal 
ballistic  parameters  do  not  vary  significantly  with  respect  to 
time  during  the  motor  operating  period,  the  difference  between  the 
results  obtained  by  these  two  methods  is  very  small.  Thus,  only 
the  results  obtained  by  the  Triple-Loop-1 teration  Algorithm  are 
presented.  For  example  2.  respective  results  obtained  by  both 
methods  are  presented. 

Examp  1 e  1 

The  grain  is  bonded  to  a  two-dimensional  hollow,  cylindrical 
case  of  star-shaped  interior.  The  motor  head  is  double-arc, 
transitionally  closed.  The  motor  tail-end  surface  is  located  on 
the  oval  portion  following  the  motor.  The  grain  combustion  is 
unrestrictive.  Calculated  results  are  shown  in  Fig.  2. 


Fig.  2.  Pressure  vs  time  for  example  1. 
Key:  cl)  Calculated  curve  (Triple-Loop- 
Iteration  Algorithm);  (2)  Test  curve. 


Examp  1 e  2 


The  grain  is  bonded  to  a  two-dimensional  hollow,  cylindrical 
case  of  wheel-shaped  interior.  Motor  head  and  tail-end  surfaces 
are  both  located  at  the  cylindrical  section  following  the  motor. 
The  grain  combustion  is  unrestrictive.  Calculated  results  are 
shown  in  Fig.  3 . 

HP 


Fig.  3.  Pressure  vs  time  for  example  2. 

Key:  1)  Calculated  curve  (Time-Space 

Alternate  Algorithm);  (2)  Calculated  curve 
(Triple-Loop- I teration  A1 gori thm) ; 

( 3 )  Test  curve . 

Example  3 

The  grain  is  of  a  single,  tube-shaped  configuration  with 
simultaneous  combustion  on  both  inner  and  outer  surfaces.  Grain 
combustion  on  both  front-  and  back-end  surfaces  are  restrictive. 
Calculated  results  are  shown  in  Fig.  4. 


Fig.  4.  Pressure  vs  time  ror  example  3. 

Key:  (1)  Calculated  curve  (Triple-Loop- 

Iteration  Algorithm);  (2)  Test  curve. 

It  can  be  seen  from  the  calculated  results  that,  when  the 
internal  ballistic  parameters  vary  considerably  with  respect  to 
time,  significant  differences  exist  between  the  results  obtained 
from  the  quasi-steady  equations  ( Tr i p 1 e-Loop- I  ter at i on  Algorithm; 
and  from  the  unsteady  equations  (Time-Space  Alternate  Algorithm). 
The  calculation  accuracy  of  the  latter  is  increased  by  about  1.6% 
compared  to  that  of  the  former. 
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STUDY  ON  ALGORITHMS  FOR  PREDICTION  OF 
SOLID  PROPELLANT  ROCKET  MOTOR  PERFORMANCE 

Wang  Xin  and  Zhang  Zhongqin 

( Beijing  Institute  of  Aeronautics  and  Astronautics ) 

Abstract 

Two  algorithms  for  prediction  of  solid  propellant  rocket  motor  perfor¬ 
mance  are  introduced.  Triple  Loop  Iteration  Algorithm  is  applied  to  solu¬ 
tion  of  quasi-steady  equationsC  5  )~(  8  ).  Its  computer  program  comprises 
a  procedure  with  three  iterative  loops.  Time-Space  Alternate  Algorithm  can 
be  applied  to  soluting  both  quasi-steady  equations  and  unsteady  equations 
(22)~(24).  If  dP/dt,  dv/dt,  dp/dt  in  equations  (22)  — ■  (24)  are  ghen,  equa¬ 
tions  (22)~(24)  will  be  reduced  to  ordinary  differential  equations.  They 
can  be  solved  by  Merson' s  method.  When  AP/Af,  Ac/ At,  Ap/Af  are  appro¬ 
ximated  to  dP/dt,  dv/dt,  dp/dt ,  this  algorithm  for  model  equation 


corresponds  to  the  finite  difference  scheme 


!“»  ut  4  “m/;  ui*i  t  t  ui~i~ Hj.i  1 


It  is  an  unconditionally  stable  scheme. 

A  general  two-dimensional  grain  calculation  method  is  introduced  for 
calculation  of  grain.  A  comprehensive  computer  program  for  above  men¬ 
tioned  methods  has  been  written  and  applied  to  calculating  the  perfor¬ 
mance  of  three  different  motors.  The  results  of  the  numerical  examples 
are  consistent  with  the  experimental  data. 
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TURBULENT  VORTICES  AND  BIONICS  IN  TURBOJET 


Ning  Huang  and  Gao  Ge 

Beijing  Institute  of  Aeronautics  and  Astronautics 
Received  2  January  1984 


The  construction  of  a  modern  turbojet  is 
rather  complex,  but  its  performance  still  has 
many  dissatisfactions,  viz.  great  pressure 
losses  along  the  flow  patch,  low  component 
efficiencies,  proneness  of  the  compressor  to 
stall'  and  surge,  frequent  occurrence  of 
overheating  and  oscillations  in  combustors, 
difficulty  in  high-altitude  reigniting,  and 
narrow  stable  operating  range.  Phenomena  of 
turbulent  vortices  are  discussed  in  this 
paper.  Newly  derived  momentum  equations  for 
turbulent  flow  are  presented.  Some  ideas 
about  morphological  simulation,  which  may 
alleviate  the  vortex  losses  and  utilize  the 
vortex  capabilities,  are  suggested  and 
thereby  produce  a  few  research  subjects  that 
are  worth  further  studying. 
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I.  Vortex  Occurs  Everywhere 

Under  various  flight  operating  conditions,  the  occurrence, 
development,  and  decay  of  turbulent  vortices  are  observed 
everywhere  along  the  flow  path  of  a  turbojet  from  its  entrance  to 
exit. 

Under  subsonic  flight  condition,  the  excursion  of  the  flow 
creates  vortices  inside  the  air  intake  of  a  turbojet. 

Under  supersonic  flight  condition,  there  are  occurrences  of 
boundary  layer  separation  or  the  laminar  flow  layers  roll  up  to 
form  turbulent  vortex  flow  field  behind  the  shock  waves  of  the  air 
intake  tunne 1 . 

In  a  compressor,  turbulent  vortices  are  created  between  the 
clearance  of  the  blade-tips  and  at  the  trailing  edge  of  the  blade. 
Between  the  stator-  and  rotor-blades,  strong  vortices  are  also 
f  ormed . 

In  an  abrupt  expansion  type  of  cylindrical  pressure  diffuser, 
stagnation  vortices  are  induced  at  the  edge  of  the  nozzle  mouth  by 
the  impinging  jet  fluid.  Fluctuations  of  the  fluid  flow  both 
upstream  and  downstream  are  dampened  because  of  the  appearance  of 
a  turbulent  flow  region  between  the  vortices  and  thereby  alleviate 
interf er ences . 

For  a  vortex  generator  or  for  injector  holes  of  a  primary 
combustion  chamber,  vortex  recirculation  region  is  created  at  the 
top  of  the  fire  box  and  thus  forms  a  heat  reservoir.  This  ensures 
an  intensified  ignition  source  and  also  shortens  the  flame. 


In  a  gas  turbine,  vortex  trails  are  installed  on  both  the 
blade  crown  and  the  reinforced  rubber  between  blades.  The  enthalpy 
drops  can  be  increased  in  a  single-stage  type  turbine  by  con¬ 
trolling  vortices  at  the  exhaust  of  the  gas  turbine.  The  number  or 
stages  is  subsequently  reduced.  This  method  has  already  been  used. 

A  forced-draft  combustion  chamber  uses  the  vortex  recircula¬ 
tion  region  behind  a  bluff-body  to  stabilize  flame. 

To  design  turbojet  components,  we  have  to  study  the  charac¬ 
teristics  of  turbulent  vortices  and  their  control  methods.  While 
vortex  turbulent  flow  consumes  energy,  loses  pressure,  reduces 
efficiency,  and  causes  malfunction,  vortices  do  possess  special 
functions  in  certain  situations. 

II.  Vortex  Characteristics 

Turbulent  flow  can  be  regarded  as  a  group  of  small-scale 

vortices.  The  vortex  characteristics  play  an  important  role  on 

« 

turbulent  flow  characteristics.  The  turbulent  flow  is  always 
present  whenever  there  are  vortices  in  the  flow  field.  The  types 
of  vortex  can  be  classified  as  vortex  filament,  vortex  sheet, 
vortex  tube,  vortex  cluster  and  eddy.  Vortex  filaments  form  a 
vortex  sheet  or  vortex  tube.  Cross-sections  of  vortex  ring  and 
vortex  cluster  are  circular  vortices.  Cross-sections  of  an  eddy 
are  spiral  vortices.  Vortex  filament  and  vortex  tube  cannot  have 
free  ends.  They  must  either  be  attached  to  the  interface  or 
connect  their  two  ends  in  order  to  form  a  vortex  cluster  or  an 
eddy.  A  vortex  sheet  cannot  exist  without  a  boundary.  They  must 
rol 1  up  to  form  an  eddy.  An  eddy  may  expand  to  a  trumpet  shape  or 
contract  to  a  cone  shape,  and  eventual ly  breaks  up  into  turbulent 
vortex  clusters  just  as  eddies  in  the  wing-tip  trails  or  cyclones 
rising  from  the  ground  surface  would  do.  The  shapes  of  vortices 


can  be  changed  by  flow  shifts  or  by  external  forces.  Their 
original  masses,  however,  remain  unchanged.  With  respect  to  its 
own  rotating  axis  or  with  respect  to  the  instantaneous  center  of 
curvature,  a  vortex  has  a  rotational  angular  velocity  co.  and  a 
circulation  T.  The  circulation  of  a  unit  cross-sectional  area  is 
called  "  vor  t  i  c  i  ty "  Q.  =  2<u.  If  the  circulation  remains  constant,  the 
product  of  the  eddy  cross-sectional  area  A  and  strength  of  vortex 
Q.  is  a  constant.  For  a  trumpet-shaped  eddy,  the  end  cross- 

sectional  area  diverges.  Thus,  cu  becomes  too  small  and  then  the 
eddy  breaks  up.  On  the  other  hand,  for  a  cone-shaped  eddy,  the 
cross-sectional  area  of  the  sharp-end  converges.  Thus,  cu  becomes 
too  large.  Since  the  micro-scale  gas  pockets  are  being  thrown  out 
by  centrifugal  force,  the  eddy  collapses  and  disperses  into  a 
turbulent  flow,  thus  uo  has  upper  and  lower  limits. 

The  circulation  of  a  circular  vortex  is  The  smaller 

the  value  of  radius  r  is,  the  faster  the  rotation  is.  An  increase 
in  strengths  of  a  circular  vortex  would  result  in  the  formation  of 
centripetal ly  converging  spiral  vortices.  The  combined  flow 
velocity  can  be  increased  from  subsonic  to  supersonic.  When  the 
limiting  flow  velocity  is  reached,  a  vacuum  vortex  will  appear:  An 
actual  vortex  core  provides  low  pressure  suction,  thus  creates 
centripetal ly  converging  flow  and  thereby  forms  a  three-dimen¬ 
sional  flow  along  its  axis.  The  vortex  core  can  be  represented  by 
a  rotating  "solid  cylinder  core"  which  induces  the  surrounding  gas 
to  follow  the  rotation.  The  tangential  velocity  at  the  inner 
circle  is  higher  than  that  at  the  outer  one.  This  is  called  an 
active  vortex.  A  mutual  induction  effect  is  observed  between 
vortices  in  this  flow  field.  Therefore,  vortex  rings  have  counter¬ 
flow,  forward-moving  velocities. 
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Shear  layers  are  generated  behind  the  trailing  eage  or  an 
obstructing  bluff-body  due  to  sudden  flow  velocity  change,  and  a 
two-dimensional  vortex  or  a  three-dimensional  votex  ring  is 
formed.  Affected  by  the  surrounding  gas  flow,  the  tangential 
velocities  inside  the  vortex  is  lower  than  that  outside  the 
vortex.  Therefore,  this  is  called  passive  vortex.  The  stability  or 
the  active  and  passive  vortices  are  different  til. 

Observations  and  experimental  studies  in  recent  years  reveal 
t  n  a  t  satisfactory  explanations,  from  the  standpoint  or  mechanism, 
tor  many  characteristics  of  vortices  are  still  waiting  to  be 
discovered.  For  example,  why  do  tornadoes  and  eddies  at  aircraft 
wing -tips  persist  for  a  long  time  before  they  decay;  why  does  the 
top  or  a  typhoon  diverge  in  reverse  rotation;  why  do  "negative 
viscosity"  regions  appear  in  the  back  section  of  the  tail  vortices 
of  a  nons t r earn  1  i ned  object  C2];  although  the  Ranque-H i 1 sch  effect 
has  already  been  applied  to  the  cold  and  hot  sublayers  of  gas 
flows,  the  explanations  of  their  mechanisms  are  still  inconclusive 
C31.  Further  understanding  of  the  fundamentals  of  vortex  charac¬ 
teristics,  therefore,  appears  to  be  a  worthwhile  research 
d i rection, 

III.  Stability  of  Vortex 

Inside  a  vortex,  flows  are  moving  in  separate  layers. 
Diffusive  exchanges  of  momentum,  energy,  and  mass  exist  when 
cross-over  flow  occurs.  Moreover,  there  is  also  viscous  dissipa¬ 
tion  and  centrifugal  force  flow  field.  Vortices  are  pressed  and 
elongated  by  the  surrounding  main  flow  or  by  the  restriction  of 
the  solid  wall  boundary  Between  the  vortex  circumferences  of  the 
recirculation  flow  and  the  co-current  boundary  layer  streamlines, 
a  narrow,  rhombus-shaped  vortex  tail  turbulent  region  is  formed 
'Fig.  2).  This  vortex  tail  turbulent  region  is  densely  populated 
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with  small  vortex  clusters  of  high  vorticity  and  is  very  unstable. 

Local  angular  velocity  63  at  the  inflection  point  of  the  long  tail 

vortex  circumference  is  very  high,  which  often  throws  gas  pockets 

in  the  recirculation  region  into  the  vortex  tail  turbulent  region 

causing  it  to  lose  stability.  If,  at  the  same  time,  there  is 

external  interference  caused  by  strong  fluid  pressure,  the  vortex 

tail  will  pull  the  vortices  loose,  which  subsequently  collapse  to 

form  a  downstream  turbulent  flow  field.  These  are  the  reasons  for 

instabilities  in  the  vortex  recirculation  region  trailing  a 

regular  vortex  generator  or  a  V-type  stabilizer.  The  authors 

believe  that  "the  flame  is  stable  only  when  the  vortex  is  stable.' 

The  stability  of  vortex  is  greatly  affected  by  the  magnitude  of 

hO 

local,  instantaneous  angular  velocity  of  vortex  and  the  length  of 
the  vortex  tail.  From  reference  Cl],  the  criterion  of  upper  and 
lower  limits  for  a  passive  vortex  is  given  by 

Jl-|+-T^-|#l<£D<|ifl  +  4-CD'  ll) 

where  {jr  gradient  of  vortex  tangential  velocity  along 
the  radius  of  curvature  r; 


frequency  of  external  interference  wave! 


co  local  angular  velocity  of  vortex; 


g  acceleration  of  a  viscous  dissipation  field!  and 


3  ratio  of  density  gradient  along  r  to  the  density 
without  interference 
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Reference  [1]  also  gives  a  criterion  of  stability  for  an 
active  vortex  and  a  spiral  vortex.  The  conditions  for  intensifying 
vortex  stability  can  be  concluded  as:  keeping  the  convergence  of 
vortex  tail  as  short  as  possible;  keeping  vortex  local  U3  within 
the  upper  and  lower  limits;  attaching  itself  to  a  solid  boundary 
to  become  a  resident  vortex;  having  the  capability  of  self-control 
such  that,  after  interference  takes  place,  the  vortex  can  be 
recovered  back  to  its  original  geometrical  shape  without  being 
detached  or  being  broken  up. 

IV.  Turbulent  Vortex  Flow  Field 

In  the  turbulent  vortex  flow  field,  the  rolling  fluctuations 
of  large  and  small  vortex  clusters  are  of  an  irregular, 
randomly-changing  process.  The  Kolmogoroff  theory  states  that 
large  vortex  clusters  disperse  and  break  up  into  small  vortex 
clusters.  Turbulent  kinetic  energy  is  gradually  transferred  from 
large  vortex  clusters  to  small  vortex  clusters.  Finally,  due  to 
viscosity,  the  turbulent  kinetic  energy  dissipates  and  results  in 
an  energy  spectrum.  Based  on  this  theory,  turbulence  model  with 
turbulent  kinetic  energy  k  and  the  turbulent  energy  dissipation  £ 
were  developed.  The  results  of  calculations  for  the  tail  flow  or  a 
bluff-body  using  the  k-£  model  show  that  the  vorticity  Q  in  the 
vortex  tail  turbulent  flow  region  is  greatly  weakened  iFig.  1). 
However,  the  experiments  in  reference  L21  verify  that  the 
vorticity  Q  in  this  region  is  very  strong  (Fig.  2). 
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Fig.  1.  Weakening  vorticity 
in  the  vortex  tail  region 
(results  of  calculated  by 
k-£  mode  1  )  . 

Key:  (l)  Boundary  layer 

streamlines;  (2)  Recircula¬ 
tion  flow  circumferences. 

In  regions  of  high  vorticity  gradient  or  pressure  sublayers, 
the  dissipation  of  local  turbulent  energy  is  reduced,  or  even 
inversed  into  macroscopic  kinetic  energ'  This  is  called  "energy- 
inversion"  or  "negative  viscosity"  phenomenon.  The  isotropic 
turbulence  model  and  statistical  methods  cannot  explain  the 
" i nterm i ttency"  and  the  seem i ng 1 y- i r r e gu 1 ar  yet  "orderly-coherent 
phenomenon"  in  a  turbulent  vortex  flow.  This  implies  that,  under 
certain  circumstances,  the  appearance  of  turbulent  vortex  clusters 
is  not  only  a  pure  transport  phenomenon  but  also  possesses 
"Soliton"  characteristics  C4,  5].  Gao  Ge  has  suggested  that  the 
canonical  equations  of  a  turbulent  vortex  flow  field  should 
include  the  "dispersion  effects"  and  combined  equations  of  the 
K-dV  equations  with  third-order  partial  derivatives  and  the 
Burgers  equations: 

+  o  ,,, 


Fig.  2.  Intensifying 
vorticity  in  the  vortex 
tail  region  (results  of 
experiments)  [23. 


where  v  is  kinematic  viscosity  coefficient  and  rj  is  dispersion 
coefficient.  The  "non- i sot r op i c" ,  "orderly-coherent  phenomenon", 
"intermittency",  and  "energy-inversion"  problems  of  a  turbulent 
flow  can  be  explained  after  an  appropriate  and  quantitative 


construction  of  the  last  two  terms  in  equation  (.Z)  have  been 
performed.  Gao  Ge  has  investigated  the  relationships  between  the 
viscous  dissipation  and  dispersion  effects.  Momentum  equations  ror 
turbulent  flow  are  derived  as: 
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where  qT  are  the  laminar  kinematic  viscosity,  turbulent 


kinematic  viscosity,  and  turbulent  dispersion  coefficient 

ax/ 


_  au:  aui 

respectively,  Vorticity  Cor  degree  of  circulation)  Q - 


The  new  momentum  equations  of  a  turbulent  flow  are  composed 
jointly  of  the  above  equations  (3)  and  (4).  Equation  c3)  is 
employed  for  low  vorticity  cases,  while  equation  (4)  is  employed 
for  high  vorticity  cases.  Details  for  determinations  of  turbulent 
kinematic  viscosity  vT  and  turbulent  dispersion  coefficient  are  the 
subject  of  another  paper.  When  the  time-averaged  vorticity  becomes 
negligibly  small,  equation  (3)  is  reduced  to  the  Boussinesq 
viscous  turbulent  vortex  equation;  when  vr  is  very  small,  the 
equation  is  the  same  as  the  laminar  Nav i er -Stokes  equation.  Manv 
turbulence  models  can  be  included  in  the  scope  of  momentum 
equat ion  ( 3 )  . 


The  orderly-coherent  phenomena  at  the  boundary  layer  and 
turbulent  vortex  flow  problems  are  successfully  solved  by  using 
the  momentum  equation  (4).  This  equation  has  been  applied  to  the 
numerical  solution  for  "the  ene r gy- i nve r s i on  phenomena  in  the 
turbulent  tail  vortex  flow  field  of  a  bluff-body".  The  calculated 
results  are  consistent  with  the  experimental  data  presented  in 
reference  C21.  Calculation  is  also  in  good  agreement  with 
experiments  of  "the  orderly-coherent  phenomena  of  a  horseshoe¬ 
shaped  pattern  in  turbulent  boundary  layers".  Reasonable  explana¬ 
tion  can  be  provided  for  many  problems  such  as  the  non-isotropic 
turbulence,  the  Ranque-Hi 1 sch  effect  of  cold  and  hot  layer 
separation,  and  the  formation  of  vortex  series  from  the  boundary 
layer  separation  on  a  wing  model,  etc. 

It  is  worthwhile  to  conduct  further  theoretical  and 
experimental  studies  on  the  dissipation  structure  and  dispersion 
phenomena  of  a  turbulent  vortex  flow  field.  For  engineering 
applications,  the  above  discussions  can  provide  design  methods  tor 
the  prevention  of  boundary  layer  separation  and  thereby  improve 
the  performance  and  efficiency  of  turbojet  components,  or  reduce 
the  drag  on  aircraft. 

V.  Turbulent  Vortices  and  Bionics 

The  interactions  between  living  things  and  land  shapes  versus 
fluid  in  the  natural  environment  take  various  forms.  As  a  fish 
swims  in  the  water,  the  vortex  trail  created  by  its  tail  fin 
movements  is  almost  indiscernible.  This  is  probably  due  to  the 
absence  of  boundary  layer  separation  or  lack  of  burden  resulted 
from  tail  vortex  drag.  Aquatic  plants  grow  densely  in  turbulent 
streams.  They  grow  in  long,  thin  bundles.  In  order  to  reduce  the 
impact  of  water  flow,  they  sway  with  the  current  and  thereby  are 
able  to  keep  from  being  uprooted  and  washed  away.  With  wings  wide 
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stretched,  an  eagle  can  glide  in  the  air  for  a  long  time.  Its 
wing-tip  feathers  open  up  like  fingers  which  weaken  wing-tip 
vortices  and  thereby  reduce  drag.  Morphological  simulation  of  this 
feature  has  been  applied  to  make  wing-tip  mounted  wingiets  on 
aircrafts.  In  the  desert,  sand  dunes  are  created  by  gusty  winds 
and  their  shapes  resemble  that  of  a  slipper's  toe.  The  front  slope 
of  a  sand  dune  is  smoothed  and  streamlined  by  the  wind,  with  its 
two  arms  extending  backward  and  touching  each  other.  Horseshoe¬ 
shaped  eddies  are  formed  behind  the  sand  dune.  This  feature 
matches  the  stability  conditions  of  a  vortex;  and  therefore,  it 
provides  the  sand  dune  with  strong  resistance  to  interference. 

The  natural  morphological  features  are  simulated  in  order  to 
improve  turbojet  components.  Some  ideas  are  first  presented  as 
foil ows : 

1)  In  order  to  simulate  the  motion  of  a  fish  mouth,  the  air 
intake  and  the  pressure  diffuser  throat  are  manufactured  such  that 
they  possess  automatic  expansion  and  contraction  type  of  struc¬ 
tures.  With  the  aid  of  signal  feedback  from  micro-sensors,  the 
peristaltic,  swallowing  capability  of  the  air  intake  tunnel  can  be 
automatically  controlled  and  thereby  prevent  surge. 

2)  In  order  to  simulate  the  swaying  motion  of  the  aquatic 
plants,  air  intake  tunnel  with  soft,  elastic  walls  may  be  manu¬ 
factured.  The  wall  surface  is  spray-coated  with  a  film  of  poly- 
tetraf luoroethylene  vinyl  coating  which  can  reduce  frictional 
force  and  deaden  pressure  fluctuations. 


40 


3)  "The  orderly-coherent  structure"  of  a  turbulent  boundary 
layer  can  be  utilized  to  design  compressor  blade  configurations. 
The  range  of  stable  compressor  operation  and  its  erficiency  can  be 
increased  by  delaying  or  preventing  the  boundary  layer  from 
separating  and  vortices  from  detaching. 

4)  In  order  to  simulate  an  eagle's  wings,  toughened  bladelets 
are  mounted  on  the  blade-tip  to  reduce  vortex  drag  losses. 

5)  "Elastic  fish  tails"  are  arranged,  along  the  height  of 
blade,  on  the  trailing  edge  of  compressor  blades  to  control 
vortex  detachment  and  reduce  tail  vortices. 

6)  A  series  of  "sand  dunes"  placed  in  a  circular  arrangement 
can  be  utilized  to  replace  vortex  recirculation  region  created  by 
vortex  generator  .in  the  head  of  a  primary  combustion  chamDer.  Such 
a  replacement  can  increase  stability  and  reduce  drag. 

7)  In  a  forced-draft  combustion  chamber,  the  phenomena  of 
hor sehsoe-shaped  eddies  being  formed  behind  a  "sand  dune"  can  be 
simulated.  This  simulation  can  result  in  increased  combustion 
stability,  oscillation  alleviation,  drag  reduction,  and  increased 
combustion  efficiency. 


References 


1.  Ge,  Gao.  On  the  Research  of  Local  Stability  of  Vortices. 
"Journal  of  Engineering  Thermophysics,”  Vol.  2,  No.  4,  p.  394, 
1981.  6th  Symposium  ISABE,  June.  1983.  Paris. 

2.  Fu.1ii,  S..  Gomi.  M.  and  Egucei.  K.  Cold  Flow  Tests  of  a 
Bluff -Body  Flame  Stabilizer.  "J.  Fluids  Engineering,"  Vol.  100, 
Sept.  1978,  p.  323. 

3.  Kurosaka,  M.  Acoustic  Streaming  in  Swirling  Flow  and 

Ranque -H i 1 sch  ( Vo r t ex -Tube )  Effect.  "J.  Fluid  Mechanics,"  Vo i . 
124,  198-2,  pp.  139-172. 

4.  Miles,  J.W.  Solitary  Waves.  "Ann.  Rev.  Fluid  Mech.  ,  "  1960,  1 

pp.  11-43. 

5.  Nishida.  Yashushi.  Experiments  on  the  Scattering  of  Ion 
Acoustic  Solitons  and  Turbulence  Formation.  "IPPJ-388,"  May,  197 


TURBULENT  VORTICES  AND  BIONICS  IN  TURBOJET 
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( Beijing  Institute  of  Aeronautics  and  Astronautics) 
Abstract 


The  construction  of  a  modern  turbojet  is  rather  complex,  but  its  perfor¬ 
mance  still  has  many  dissatisfactions,  viz.  great  pressure  losses  along  the 
flow  path,  low  component  efficiencies,  narrow  stable  operating  range, 
proneness  of  the  compressor  to  stall  and  surge,  difficulty  in  high-altitude 
relighting,  frequent  occurrence  of  overheating  and  oscillations  in  combustors. 
The  authors  consider  that  all  the  above  drawbacks  are  closely  related  to 
turbulent  vortices,  which  occur  almost  everywhere  within  and  around 
a  turbojet.  The  generation,  development  and  decay  of  vortices  have  profound 
mutual  effects  with  the  surrounding  flow.  In  aerodynamic  design  an  endea¬ 
vour  must  be  made  to  reduce  vortex  losses  and  exploit  vertex  capabilities. 
There  are  some  interesting  peculiar  phenomena  in  a  turbulent  vortex  flow, 
such  as  orderly-coherent  flow  pattern,  intermittency  and  energy-inversion 
unconformable  to  the  Kolmogoroff  theory  of  energy  spectrum.  The  tradi¬ 
tional  isotropic  turbulence  and  statistical  treatment  cannot  explain  the  above 
phenomena.  The  co-author  has  discovered  that  viscous  dissipation  and  dis¬ 
persion  effects  should  be  considered  simultaneously  and  Burger's  and  the 
K_dV  equations  have  to  be  mixed  up  to  obtain  a  canonical  equation! 
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According  to  this  concept  two  new  momentum  equations  of  turbulent  flow 
are  derived: 
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Sample  computations  of  a  tail  vortex  flow  and  an  orderly-coherent  flow 
pattern  have  been  performed  with  success.  The  authors  suggest  some  ideas 
about  morphological  simulation  which  may  alleviate  the  vortex  evils  and 
utilize  the  vortex  characteristics. 


